<html>
 <head>
  <meta charset="utf-8"/>
  <meta content="width=device-width, initial-scale=1, maximum-scale=1, user-scalable=no" name="viewport"/>
  <title>
   邓一硕: 蒙特卡洛方法与定积分计算  | 数螺 | NAUT IDEA
  </title>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap-theme.min.css" rel="stylesheet"/>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap.min.css" rel="stylesheet"/>
  <style type="text/css">
   #xmain img {
                  max-width: 100%;
                  display: block;
                  margin-top: 10px;
                  margin-bottom: 10px;
                }

                #xmain p {
                    line-height:150%;
                    font-size: 16px;
                    margin-top: 20px;
                }

                #xmain h2 {
                    font-size: 24px;
                }

                #xmain h3 {
                    font-size: 20px;
                }

                #xmain h4 {
                    font-size: 18px;
                }


                .header {
	           background-color: #0099ff;
	           color: #ffffff;
	           margin-bottom: 20px;
	        }

	        .header p {
                  margin: 0px;
                  padding: 10px 0;
                  display: inline-block;  
                  vertical-align: middle;
                  font-size: 16px;
               }

               .header a {
                 color: white;
               }

              .header img {
                 height: 25px;
              }
  </style>
  <script src="http://cdn.bootcss.com/jquery/3.0.0/jquery.min.js">
  </script>
  <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" type="text/javascript">
   MathJax.Hub.Config({
          tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]},
          TeX: {equationNumbers: {autoNumber: "AMS"}}
        });
  </script>
  <script src="http://nautstatic-10007657.file.myqcloud.com/static/css/readability.min.js" type="text/javascript">
  </script>
  <script type="text/javascript">
   $(document).ready(function() {
                 var loc = document.location;
                 var uri = {
                  spec: "http://cos.name/2010/03/monte-carlo-method-to-compute-integration/",
                  host: "http://cos.name",
                  prePath: "http://cos.name",
                  scheme: "http",
                  pathBase: "http://cos.name/"
                 };
    
                 var documentClone = document.cloneNode(true);
                 var article = new Readability(uri, documentClone).parse();
     
                 document.getElementById("xmain").innerHTML = article.content;
                });
  </script>
  <!-- 1466441954: Accept with keywords: (title(0.25):门户网站,统计学,计算,服务平台,中国,蒙特卡洛,定,方法, topn(0.5):稀疏,R语言,事件,目录,函数,模拟法,平均值,结果,均匀分布,维数,例子,随机数,文章,蒙特卡罗,数据,抽样,投点,计算,样本,问题,命令,算法,随机,区间,方法,代码,评论,版本,时候,蒙特卡洛).-->
 </head>
 <body class="single single-post postid-2010 single-format-standard sidebar" onload="">
  <div class="header">
   <div class="container">
    <div class="row">
     <div class="col-xs-6 col-sm-6 text-left">
      <a href="/databee">
       <img src="http://nautidea-10007657.cos.myqcloud.com/logo_white.png"/>
      </a>
      <a href="/databee">
       <p>
        数螺
       </p>
      </a>
     </div>
     <div class="hidden-xs col-sm-6 text-right">
      <p>
       致力于数据科学的推广和知识传播
      </p>
     </div>
    </div>
   </div>
  </div>
  <div class="container text-center">
   <h1>
    邓一硕: 蒙特卡洛方法与定积分计算
   </h1>
  </div>
  <div class="container" id="xmain">
   <div class="hfeed site" id="page">
    <header class="site-header" id="masthead" role="banner">
     <div id="cos-logo">
      <a href="http://cos.name/">
       <img src="http://cos.name/wp-content/themes/COS-kermesinus/images/headers/cos-logo.png"/>
      </a>
     </div>
     <div class="navbar" id="navbar">
      <nav class="navigation main-navigation" id="site-navigation" role="navigation">
       <h3 class="menu-toggle">
        菜单
       </h3>
       <a class="screen-reader-text skip-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#content" title="跳至内容">
        跳至内容
       </a>
       <div class="menu-%e6%88%91%e7%9a%84%e8%8f%9c%e5%8d%95-container">
        <ul class="nav-menu" id="menu-%e6%88%91%e7%9a%84%e8%8f%9c%e5%8d%95">
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-home menu-item-4746" id="menu-item-4746">
          <a href="http://cos.name">
           主页
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-has-children menu-item-8120" id="menu-item-8120">
          <a href="http://cos.name/cn">
           论坛
          </a>
          <ul class="sub-menu">
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8124" id="menu-item-8124">
            <a href="http://cos.name/cn/wp-login.php?action=register">
             论坛注册
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8125" id="menu-item-8125">
            <a href="http://cos.name/cn/wp-login.php">
             论坛登录
            </a>
           </li>
          </ul>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8110" id="menu-item-8110">
          <a href="http://cos.name/books/">
           图书资料
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8112" id="menu-item-8112">
          <a href="http://cos.name/videos">
           视频教程
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8333" id="menu-item-8333">
          <a href="http://cos.name/salon/">
           统计沙龙
          </a>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8111" id="menu-item-8111">
          <a href="http://cos.name/chinar/">
           R语言会议
          </a>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8109" id="menu-item-8109">
          <a href="http://cos.name/training/">
           讲座与培训
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8115" id="menu-item-8115">
          <a href="http://cos.name/cn/forum/comprehensive/job/">
           招聘信息
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-has-children menu-item-4780" id="menu-item-4780">
          <a href="http://cos.name/about">
           关于我们
          </a>
          <ul class="sub-menu">
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8113" id="menu-item-8113">
            <a href="http://cos.name/2008/11/how-to-work-with-cos/">
             加入我们
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8119" id="menu-item-8119">
            <a href="http://cos.name/donate/">
             赞助我们
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8114" id="menu-item-8114">
            <a href="http://cos.name/cn/forum/26">
             项目合作
            </a>
           </li>
          </ul>
         </li>
        </ul>
       </div>
       <form action="http://cos.name/" class="search-form" method="get" role="search">
        <label>
         <span class="screen-reader-text">
          搜索：
         </span>
        </label>
       </form>
      </nav>
      <!-- #site-navigation -->
     </div>
     <!-- #navbar -->
    </header>
    <!-- #masthead -->
    <div class="site-main" id="main">
     <div class="content-area" id="primary">
      <div class="site-content" id="content" role="main">
       <article class="post-2010 post type-post status-publish format-standard hentry category-optim-sim tag-monte-carlo tag-320 tag-244 tag-321" id="post-2010">
        <header class="entry-header">
         <h1 class="entry-title">
          蒙特卡洛方法与定积分计算
         </h1>
         <div class="entry-meta">
          <span class="date">
           <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/" rel="bookmark" title="链向蒙特卡洛方法与定积分计算的固定链接">
            <time class="entry-date" datetime="2010-03-08T17:21:23+00:00">
             2010/03/08
            </time>
           </a>
          </span>
          <span class="categories-links">
           <a href="http://cos.name/category/computing/optim-sim/" rel="category tag">
            优化与模拟
           </a>
          </span>
          <span class="tags-links">
           <a href="http://cos.name/tag/monte-carlo/" rel="tag">
            Monte-Carlo
           </a>
           、
           <a href="http://cos.name/tag/%e5%ae%9a%e7%a7%af%e5%88%86/" rel="tag">
            定积分
           </a>
           、
           <a href="http://cos.name/tag/%e6%a8%a1%e6%8b%9f/" rel="tag">
            模拟
           </a>
           、
           <a href="http://cos.name/tag/%e8%92%99%e7%89%b9%e5%8d%a1%e6%b4%9b/" rel="tag">
            蒙特卡洛
           </a>
          </span>
          <span class="author vcard">
           <a class="url fn n" href="http://cos.name/author/dengyishuo/" rel="author" title="查看所有由邓一硕发布的文章">
            邓一硕
           </a>
          </span>
         </div>
         <!-- .entry-meta -->
        </header>
        <!-- .entry-header -->
        <div class="entry-content">
         <p>
          本文讲述一下蒙特卡洛模拟方法与定积分计算，首先从一个题目开始：设$0\leq f(x) \leq 1$，用蒙特卡洛模拟法求定积分$J=\int_{0}^{1}f(x)dx$的值。
         </p>
         <p style="text-align: center;">
          <strong>
           随机投点法
          </strong>
         </p>
         <p>
          设$(X,Y)$服从正方形 $\{0\leq x \leq 1,0\leq y\leq 1\}$上的均匀分布，则可知 $X,Y$分别服从[0,1]上的均匀分布，且$X,Y$相互独立。记事件$A=\{Y\leq f(X)\}$，则$A$的概率为
         </p>
         <p>
          $
          <br/>
          P(A)=P(Y\leq f(X))=\int_{0}^{1}\int_{0}^{f(x)}dydx=\int_{0}^{1}f(x)dx=J
          <br/>
          $
         </p>
         <p>
          即定积分$J$的值 就是事件$A$出现的频率。同时，由伯努利大数定律，我们可以用重复试验中$A$出现的频率作为 $p$的估计值。即将$(X,Y)$看成是正方形$\{0\leq x \leq 1,0\leq y \leq 1\}$内的随机投点，用随机点落在区域${y\leq f(x)}$中的频率作为定积分的近似值。这种方法就叫随机投点法，具体做法如下：
         </p>
         <p style="text-align: center;">
          <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/m-6/" rel="attachment wp-att-2051">
           <img src="http://cos.name/wp-content/uploads/2010/03/m5.png"/>
          </a>
         </p>
         <p style="text-align: center;">
          <a href="http://cos.name/../../wp-content/uploads/2010/03/m4.png">
          </a>
          图1 随机投点法示意图
         </p>
         <p>
          1、首先产生服从 $[0,1]$上的均匀分布的$2n$个随机数（ $n $为随机投点个数，可以取很大，如 $n=10^4$ ）并将其配对。
         </p>
         <p>
          2、对这$n$对数据 $(x_i,y_i),i=1,2,…,n$ ，记录满足不等式$y_i\leq f(x_i)$的个数，这就是事件 $A$ 发生的频数$\mu_n$，由此可得事件$A$发生的频率$\frac{\mu_n}{n}$，则$J\approx \frac{\mu_n} {n}$。
         </p>
         <p>
          举一实例，譬如要计算$\int_{0}^{1}e^{-x^2/2}/\sqrt{2\pi}dx$，模拟次数$n=10^4$时，R代码如下：
         </p>
         <pre class="brush: r">n=10^4;
 x=runif(n);
 y=runif(n);
 f=function(x)
 {
 exp(-x^2/2)/sqrt(2*pi)
 }
 mu_n=sum(y&lt;f(x));
 J=mu_n/n;
 J</pre>
         <p>
          模拟次数 $n=10^5$ 时，令$n=10^5$,其余不变。
         </p>
         <p>
          定积分$\int_{0}^{1}e^{-x^2/2}/\sqrt{2\pi}dx$的精确值和模拟值如下：
         </p>
         <p>
          表1
         </p>
         <table border="1" cellpadding="0" cellspacing="0">
          <tbody>
           <tr>
            <td valign="top">
             精确值
            </td>
            <td valign="top">
             $10^3$
            </td>
            <td valign="top">
             $10^4$
            </td>
            <td valign="top">
             $10^5$
            </td>
            <td valign="top">
             $10^6$
            </td>
            <td valign="top">
             $10^7$
            </td>
           </tr>
           <tr>
            <td valign="top">
             0.3413447
            </td>
            <td valign="top">
             0.342
            </td>
            <td valign="top">
             0.344
            </td>
            <td valign="top">
             0.34187
            </td>
            <td valign="top">
             0.341539
            </td>
            <td valign="top">
             0.341302
            </td>
           </tr>
          </tbody>
         </table>
         <pre class="brush: r">注：精确值用integrate(f,0,1)求得</pre>
         <p>
          扩展
         </p>
         <p>
          如果你很细心，你会发现这个方法目前只适用于积分区间$[0,1]$ ，且积分函数 $f(x)$ 在区间$[0,1]$上的取值也位于 $[0,1]$ 内的情况。那么，对于一般区间 $[a,b]$ 上的定积分$J’=\int_{a}^{b}g(x)dx$ 呢？一个很明显的思路，如果我们可以将 $J’$ 与$J$ 建立代数关系就可以了。
         </p>
         <p>
          首先，做线性变换，令 $y=(x-a)/(b-a)$ ，此时，
         </p>
         <p>
          $x=(b-a)y+a$, $J’=(b-a)\int_{0}^{1}g[(b-a)y+a]dy$。
         </p>
         <p>
          进一步如果在区间$[a,b]$上有$c\leq g(x) \leq d $ ，令
         </p>
         <p>
          $f(y)=\frac{1}{d-c}{g(x)-c}=\frac{1}{d-c}{g[a+(b-a)y]-c} $，
         </p>
         <p>
          则$0\leq f(y) \leq 1$。此时，可以得到$J’=\int_{a}^{b}g(x)dx=S_0J+c(b-a) $。
         </p>
         <p>
          其中,$S_0=(b-a)(d-c) $, $J=\int_{0}^{1}f(y)dy $。
         </p>
         <p>
          这说明，用随机投点法计算定积分方法具有普遍意义。
         </p>
         <p>
          举一个实例，求定积分 $J’=\int_{2}^{5}e^{-x^2/2}/\sqrt{2\pi}dx $。
         </p>
         <p>
          显然$a=2,b=5$，$c=min\{g(x)|2\leq x \leq 5\},d=max\{g(x)|2\leq x \leq 5\}$，由于$g(x)=e^{-x^2/2}/\sqrt{2\pi}$在 $[2,5]$上时单调减函数，所以$c=g(5),d=g(2)$，$S_0=(b-a)(d-c)$。R中代码为
         </p>
         <pre class="brush: r">a=2;
 b=5;
 g=function(x)
 {
 exp(-x^2/2)/sqrt(2*pi);
 }
 f=function(y)
 {
 (g(a+(b-a)*y)-c)/(d-c);
 }
 c=g(5);d=g(2);s_0=(b-a)*(d-c);
 n=10^4;
 x=runif(n);y=runif(n);
 mu_n=sum(y&lt;=f(x));
 J=mu_n/n;
 J_0=s_0*J+c*(b-a);</pre>
         <p>
          定积分 $J’=\int_{2}^{5}e^{-x^2/2}/\sqrt{2\pi}dx$ 的精确值和模拟值如下：
         </p>
         <p>
          表2
         </p>
         <table border="1" cellpadding="0" cellspacing="0">
          <tbody>
           <tr>
            <td valign="top">
             真实值
            </td>
            <td valign="top">
             $10^3$
            </td>
            <td valign="top">
             $10^4$
            </td>
            <td valign="top">
             $10^5$
            </td>
            <td valign="top">
             $10^6$
            </td>
            <td valign="top">
             $10^7$
            </td>
           </tr>
           <tr>
            <td valign="top">
             0.02274985
            </td>
            <td valign="top">
             0.02332792
            </td>
            <td valign="top">
             0.02311736
            </td>
            <td valign="top">
             0.02262659
            </td>
            <td valign="top">
             0.02284152
            </td>
            <td valign="top">
             0.02278524
            </td>
           </tr>
          </tbody>
         </table>
         <pre class="brush: r">注：精确值用integrate(g,2,5)求得)</pre>
         <p>
          平均值法
         </p>
         <p>
          蒙特卡洛模拟法计算定积分时还有另一种方法，叫平均值法。这个原理也很简单：设随机变量 $X$ 服从$[0,1]$上的均匀分布，则$Y=f(X)$的数学期望为
         </p>
         <p>
          $E(f(X))=\int_{0}^{1}f(x)dx=J $
         </p>
         <p>
          所以估计$J$的值就是估计$f(X)$的数学期望值。由辛钦大数定律，可以用$f(X)$的观察值的均值取估计$f(X)$的数学期望。具体做法：
         </p>
         <p>
          先用计算机产生$n$个服从$[0,1]$上均匀分布的随机数：$x_i,i=1,2,…,n$。
         </p>
         <p>
          对每一个$x_i$ ，计算$f(x_i)$。
         </p>
         <p>
          计算$J\approx \frac{1}{n}\sum_{i=1}^{n}f(x_i)$。
         </p>
         <p>
          譬如，计算 $J=\int_{0}^{1}e^{-x^2/2}/\sqrt{2\pi}dx $，R中的代码为
         </p>
         <pre class="brush: r">n=10^4;
 x=runif(n);
 f=function(x)
 {
 exp(-x^2/2)/sqrt(2*pi)
 }
 J=mean(f(x));</pre>
         <p>
          其精确值和模拟值如下：
          <br/>
          表3
         </p>
         <table border="1" cellpadding="0" cellspacing="0">
          <tbody>
           <tr>
            <td valign="top">
             真实值
            </td>
            <td valign="top">
             $10^3$
            </td>
            <td valign="top">
             $10^4$
            </td>
            <td valign="top">
             $10^5$
            </td>
            <td valign="top">
             $10^6$
            </td>
            <td valign="top">
             $10^7$
            </td>
           </tr>
           <tr>
            <td valign="top">
             0.3413447
            </td>
            <td valign="top">
             0.3405831
            </td>
            <td valign="top">
             0.3410739
            </td>
            <td valign="top">
             0.3414443
            </td>
            <td valign="top">
             0.3414066
            </td>
            <td valign="top">
             0.3413366
            </td>
           </tr>
          </tbody>
         </table>
         <p>
          平均值法与随机投点法类似可以进行扩展，这里不再赘述。
         </p>
         <p style="text-align: center;">
          <strong>
           结论
          </strong>
         </p>
         <p>
          用蒙特卡洛模拟法计算定积分具有普遍意义。蒙特卡洛模拟方法为我们提供了一个看待世界的新思路，即一个不具随机性的事件可以通过一定的方法用随机事件来模拟或逼近。
         </p>
         <div class="wumii-hook">
          <br/>
          <br/>
         </div>
        </div>
        <!-- .entry-content -->
        <footer class="entry-meta">
         <div class="author-info">
          <div class="author-avatar">
           <img src="http://sdn.geekzu.org/avatar/7969a4509f4093e8637e6bf70f3f7fc1?s=74&amp;d=monsterid&amp;r=g"/>
          </div>
          <!-- .author-avatar -->
          <div class="author-description">
           <h2 class="author-title">
            关于邓一硕
           </h2>
           <p class="author-bio">
            专注于R语言在金融投资分析和计量经济学中的应用；《R in a Nut Shell》、《R Graphics Cookbook》、《Introductory Statistics with R》等书的译者。
            <a class="author-link" href="http://cos.name/author/dengyishuo/" rel="author">
             查看所有由邓一硕发表的文章
             <span class="meta-nav">
              →
             </span>
            </a>
           </p>
          </div>
          <!-- .author-description -->
         </div>
         <!-- .author-info -->
        </footer>
        <!-- .entry-meta -->
       </article>
       <!-- #post -->
       <nav class="navigation post-navigation" role="navigation">
        <h1 class="screen-reader-text">
         文章导航
        </h1>
        <div class="nav-links">
         <a href="http://cos.name/2010/01/generalization-to-buffons-needle/" rel="prev">
          <span class="meta-nav">
           ←
          </span>
          蒲丰投针问题的推广
         </a>
         <a href="http://cos.name/2010/04/kernel-density-estimation-with-unbounded-region/" rel="next">
          有边界区间上的核密度估计
          <span class="meta-nav">
           →
          </span>
         </a>
        </div>
        <!-- .nav-links -->
       </nav>
       <!-- .navigation -->
       <div class="comments-area" id="comments">
        <h2 class="comments-title">
         《
         <span>
          蒙特卡洛方法与定积分计算
         </span>
         》有38个想法
        </h2>
        <ol class="comment-list">
         <li class="comment even thread-even depth-1" id="comment-1239">
          <article class="comment-body" id="div-comment-1239">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/34e0143fb09d6f1a82876586307be50c?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              huang shuai
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1239">
              <time datetime="2010-03-10T15:34:51+00:00">
               2010/03/10 15:34
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             沙发~~
             <br/>
             蒙特卡洛法是现代统计计算的基石，是贝叶斯获第二春的精髓所在，只是可惜，遇到高维向量，蒙特卡洛法就死翘翘了。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给huang shuai" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1239#respond" onclick='return addComment.moveForm( "div-comment-1239", "1239", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-odd thread-alt depth-1 parent" id="comment-1242">
          <article class="comment-body" id="div-comment-1242">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/5c85d270961fb44f3b10bd684f0c60ee?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              焦静
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1242">
              <time datetime="2010-03-10T19:06:27+00:00">
               2010/03/10 19:06
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             以上的程序是二维的，如果是三维的，比如这样一个积分，相当于求一个曲面y+z=1在x,y,z都在[0,1]间上积分，那么可以明显肉眼测得，这个积分值为0.5，如果用蒙特卡洛的思路，那么：
            </p>
            <pre class="brush: r">n1 &lt; - scan()
n &lt;- as.numeric(n1)
x &lt;- runif(n)
y &lt;- runif(n)
z &lt;- runif(n)
z1 &lt;- function(y) {
    1 - y
}
num &lt;- sum(z &lt; z1(y))
J &lt;- num/n
J 
</pre>
            <p>
             n1是自己输入的一个很大的数，
             <br/>
             这样的J也是在0.5附近的，我模拟的结果是0.500552
             <br/>
             其中的n1设成1000000，所以，我觉得有些情况，这种
             <br/>
             方法如果可以降维的话，也可以解决一些高维度问题吧！
             <br/>
             PS：自己乱想的，献丑了！
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给焦静" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1242#respond" onclick='return addComment.moveForm( "div-comment-1242", "1242", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment byuser comment-author-yihui even depth-2 parent" id="comment-1247">
            <article class="comment-body" id="div-comment-1247">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                <a class="url" href="http://yihui.name" rel="external nofollow">
                 谢益辉
                </a>
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1247">
                <time datetime="2010-03-11T05:19:05+00:00">
                 2010/03/11 05:19
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               高维的情况下马上就遇到维数灾难了，即使样本点非常多，在高维空间下也是很稀疏的，所以蒙特卡洛不会太精确。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给谢益辉" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1247#respond" onclick='return addComment.moveForm( "div-comment-1247", "1247", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
            <ol class="children">
             <li class="comment odd alt depth-3" id="comment-1251">
              <article class="comment-body" id="div-comment-1251">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/5c85d270961fb44f3b10bd684f0c60ee?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  焦静
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1251">
                  <time datetime="2010-03-11T11:16:55+00:00">
                   2010/03/11 11:16
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 哦，好像还是不是很了解，我再研究研究吧，高维用降维思路怎么就不可以了啊，哎，本人愚笨，要不小谢给我举个例子吧！另外，我最近才算是踏踏实实的准备从新扫盲学习R，但执行你上面的程序时出现了问题，百思不得其解啊！特来请教这个低级问题：我用2.8.1版本的R去加载animation版本1.0-4后，我就用了一个命令update.packages()，这个命令不是能将所有安装包都更新么，但是后面继续执行：
                 <br/>
                 library(animation)
                 <br/>
                 ani.options(interval = 0.2, nmax = 100)
                 <br/>
                 MC.hitormiss()
                 <br/>
                 后，仍然报错为：错误: 没有”MC.hitormiss”这个函数
                 <br/>
                 我直接用update.packages(“animation”)的命令也总是报错为：
                </p>
                <p>
                 Warning message:
                 <br/>
                 In list.files(lib) : list.files: ‘animation’目录不可读
                </p>
                <p>
                 而我同学直接安装后，就可以直接执行以上命令，他的R是2.10.1的，
                 <br/>
                 所以我们推测可以是我的R版本太低的缘故，是这样么？如果是，那么
                 <br/>
                 我还想请教小谢，怎么可以在不改变R版本的情况下，将library里面
                 <br/>
                 特定的包更新啊？
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-yihui even depth-3" id="comment-1256">
              <article class="comment-body" id="div-comment-1256">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://yihui.name" rel="external nofollow">
                   谢益辉
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1256">
                  <time datetime="2010-03-11T12:48:28+00:00">
                   2010/03/11 12:48
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 所谓维数灾难一般也就是指高维空间下数据的稀疏性。这个很容易理解，比如有10000个样本点均匀分布在样本空间里，如果数据的维数是2，那么正方形的每条边平均有100个样本点；如果数据维数是4，那么一个4维“正方体”每条边上平均只能分布10个数据点；如果维数更高，那么每一维都只能分布到更少的数据点。
                </p>
                <p>
                 在蒙特卡洛积分这里，如果是高维积分，那么抽取随机数的时候就很难照顾到每个变量，也就是说每个变量根本抽不到几个数字样本量就已经很大了，这样由于随机性很差（试想rnorm(2)能算是正态分布的代表么？），积分的结果可能就不靠谱。不知道我这种解释是不是对的。
                </p>
                <p>
                 用程序来说，计算一个n重积分[latex]\int_0^1\int_0^1\cdots\int_0^1x_1x_2\cdots x_ndx_1dx_2\cdots dx_n[/latex]，这个积分的结果我们知道，就是[latex]1/2^n[/latex]，用蒙特卡洛积分算算模拟的值，并和真实值比较：
                </p>
                <pre class="brush: r"># 程序写得稍微有点晦涩，大概意思：
# 10个均匀分布随机数相乘（x1x2...xn），看看比1个均匀分布
# （y）小的概率，考虑到数值精度问题，把乘法改成了exp(log(加法))
# ndim：数据维数；nrep：抽样次数
intProd = function(ndim, nrep = 1e+05) {
    c(MC.value = mean(runif(nrep, 0, 1) &lt; = exp(replicate(nrep, 
        sum(log(runif(ndim)))))), TRUE.value = 1/2^ndim)
}
set.seed(123)
res1 = replicate(10, intProd(15))
res2 = replicate(10, intProd(1))
# 模拟值/真实值
&gt; res1[1, ]/res1[2, ]
 [1] 1.96608 0.65536 0.98304 1.31072 1.63840
 [6] 0.65536 1.63840 1.63840 1.96608 0.32768
&gt; res2[1, ]/res2[2, ]
 [1] 1.00626 1.00260 1.00078 0.99894 0.99892
 [6] 1.00616 1.00214 0.99872 1.00050 1.00618
</pre>
                <p>
                 可以看出，样本量固定为十万，当维数为15的时候，模拟结果和真实结果的倍数波动相对比较大，维数为1的时候，结果波动相对较小。
                </p>
                <p>
                 至于animation包的问题，你最好是装最新版的R。不然的话就在加载animation包之前运行update.packages()，否则包的目录被R进程占用了就不能再安装写入了。
                </p>
                <p>
                 如果用Windows，装R的时候我有个习惯，就是把安装过程中的安装目录的版本号去掉，直接装到C:\Program Files\R\目录下，而不带R-2.8.1之类的子目录。这样每次都安装到同一个目录下，卸载R的时候不会卸载以前安装的附加包；在重新安装R之后，只需要update.packages()就可以更新所有可更新的包了。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment odd alt depth-3" id="comment-1263">
              <article class="comment-body" id="div-comment-1263">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/e787625af13288963c5a63bd2a4e5e82?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://gaotao.name" rel="external nofollow">
                   gaotao
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1263">
                  <time datetime="2010-03-12T00:35:57+00:00">
                   2010/03/12 00:35
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 冒个泡，同意谢老大的观点。因为随机投点法比较粗糙，而且用的随机数应该是伪随机的，遇到高维时，单纯的加大投点的随机数，意义不是很大，当然还因为在高维下数据太稀疏的缘故，所以随机投点法虽常用但是对于需要比较精确一点的解时误差还是挺大的。
                 <br/>
                 而对于平均值法，现在改进的方法倒是还蛮多的。在平均值法中的均方误差等于E|e_N|^2=Var(f)/N，Var(f)为f(X)的方差，所以改进时需要增大N和减小Var(f)，因此抽样方法的选择就是一个很重要的技巧，这里就不再过多探讨，可以参阅相关计算方法的书。用这个方法的话，根据上面的公式可知MC收敛阶不依赖维数（半阶收敛性），因此用较大型的计算机便可以得到较精确的解，而不用随机投点法了。
                 <br/>
                 但是当有更高效率和精度的确定性算法时，还是不会用MC法的，他是往往没有更好的算法时才采用，比如上千上万维的积分…
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment even depth-3" id="comment-1305">
              <article class="comment-body" id="div-comment-1305">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/6914198801aec79d4aae8c3d855869d3?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  胡荣兴
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1305">
                  <time datetime="2010-03-30T20:55:25+00:00">
                   2010/03/30 20:55
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 很久之前用过蒙特卡罗方法，现在都生疏了。蒙特卡罗的好处就是计算简化，不用去考虑那么复杂的积分或微分甚至偏微分。
                 <br/>
                 而对维数不敏感也是其一个优点之一，在谢兄这倒成了其缺点了。有两位经济学家用该方法解作过1900多维计算。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-yihui odd alt depth-3" id="comment-1306">
              <article class="comment-body" id="div-comment-1306">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://yihui.name" rel="external nofollow">
                   谢益辉
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1306">
                  <time datetime="2010-03-30T22:52:47+00:00">
                   2010/03/30 22:52
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 那要看怎么理解“对维数不敏感”了，如果只是“不限制维数”的话，那这算不上优点，因为高维情况下需要的样本点会急剧增加，否则积分可能不精确（如前例）。当然，蒙特卡洛方法包括思想类似的一大类方法，也许存在某些方法能克服维数问题，我不了解。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
            </ol>
            <!-- .children -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1" id="comment-1243">
          <article class="comment-body" id="div-comment-1243">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/5c85d270961fb44f3b10bd684f0c60ee?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              焦静
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1243">
              <time datetime="2010-03-10T19:08:16+00:00">
               2010/03/10 19:08
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             咦，我先试试看这次的内容可否提交，怎么我刚才提交的东西没有了呢？
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给焦静" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1243#respond" onclick='return addComment.moveForm( "div-comment-1243", "1243", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-dengyishuo bypostauthor odd alt thread-odd thread-alt depth-1 parent" id="comment-1244">
          <article class="comment-body" id="div-comment-1244">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/7969a4509f4093e8637e6bf70f3f7fc1?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              邓一硕
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1244">
              <time datetime="2010-03-10T23:24:55+00:00">
               2010/03/10 23:24
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             有的，我看得到。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给邓一硕" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1244#respond" onclick='return addComment.moveForm( "div-comment-1244", "1244", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment byuser comment-author-yihui even depth-2 parent" id="comment-1245">
            <article class="comment-body" id="div-comment-1245">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                <a class="url" href="http://yihui.name" rel="external nofollow">
                 谢益辉
                </a>
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1245">
                <time datetime="2010-03-10T23:52:11+00:00">
                 2010/03/10 23:52
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               刚才被误判为垃圾评论了。已经恢复回来。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给谢益辉" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1245#respond" onclick='return addComment.moveForm( "div-comment-1245", "1245", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
            <ol class="children">
             <li class="comment odd alt depth-3" id="comment-1262">
              <article class="comment-body" id="div-comment-1262">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/5c85d270961fb44f3b10bd684f0c60ee?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  焦静
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1262">
                  <time datetime="2010-03-11T21:03:59+00:00">
                   2010/03/11 21:03
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 多谢小谢了，这个例子举得很好，我一看就懂了你的意思了。另外，
                 <br/>
                 多维情况的那个例子中，有个命令set.seed(123)，我发现如果
                 <br/>
                 我不执行这句程序，也可以得到结果，请教了几个同学，最后我们认为
                 <br/>
                 set.seed的作用是不是可以防止重复运行intProd时候，产生相同的随
                 <br/>
                 机数？以及里面参数123是如何设定的啊？是不是需要遵从什么规则啊？
                </p>
                <p>
                 对了，我还想请教：我发现很多时候library和require使用上没什么差别，
                 <br/>
                 看了help后也发现不了什么差别，难道真的没有差别？为什么R不用一个而是
                 <br/>
                 存在两个一样的命令呢？
                 <br/>
                 PS：我初学，问题多多。劳烦小谢了。
                 <br/>
                 关于安装R包的问题，非常感谢小谢提供的意见，我就按着你说的做，太好
                 <br/>
                 了，方便了很多，谢谢。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-yixuanq even depth-3" id="comment-1272">
              <article class="comment-body" id="div-comment-1272">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/f09663fd5fc04bad5e5b09daddb86fe6?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  邱怡轩
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1272">
                  <time datetime="2010-03-14T09:07:11+00:00">
                   2010/03/14 09:07
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 <code>
                  set.seed()
                 </code>
                 不是为了防止产生相同的随机数，相反，正是为了产生相同的随机数，因为这样一来，每次运行的结果就是一样的，这使得结果可重复被验证。参数的选择很随意，比如选个幸运数字什么的。
                 <br/>
                 <code>
                  library()
                 </code>
                 和
                 <code>
                  require()
                 </code>
                 在作用上没什么大差别，但
                 <code>
                  require()
                 </code>
                 会返回一个逻辑值，加载成功就返回TRUE，否则为FALSE。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
            </ol>
            <!-- .children -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-yihui odd alt thread-even depth-1" id="comment-1246">
          <article class="comment-body" id="div-comment-1246">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              <a class="url" href="http://yihui.name" rel="external nofollow">
               谢益辉
              </a>
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1246">
              <time datetime="2010-03-11T05:12:14+00:00">
               2010/03/11 05:12
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             其实你可以做两个简单的动画图片插进来，这样会让文章更生动。你写的这两种方法我都已经在animation包中实现了：
            </p>
            <pre class="brush: r">## 需要最新版本的animation包（&gt;=1.0-10）
library(animation)
ani.options(interval = 0.2, nmax = 100)

## 随机投点
MC.hitormiss()

## 取随机数，算函数实现，算样本均值
MC.samplemean()
</pre>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给谢益辉" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1246#respond" onclick='return addComment.moveForm( "div-comment-1246", "1246", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-yihui even thread-odd thread-alt depth-1 parent" id="comment-1248">
          <article class="comment-body" id="div-comment-1248">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              <a class="url" href="http://yihui.name" rel="external nofollow">
               谢益辉
              </a>
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1248">
              <time datetime="2010-03-11T06:05:00+00:00">
               2010/03/11 06:05
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             第一个例子中精确值是用
             <code>
              pnorm(1) - pnorm(0)
             </code>
             算出来的么？如果是，可以注明一下，因为这个积分似乎没有显式表达式（正态分布的密度嘛）。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给谢益辉" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1248#respond" onclick='return addComment.moveForm( "div-comment-1248", "1248", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment byuser comment-author-dengyishuo bypostauthor odd alt depth-2 parent" id="comment-1250">
            <article class="comment-body" id="div-comment-1250">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/7969a4509f4093e8637e6bf70f3f7fc1?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                邓一硕
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1250">
                <time datetime="2010-03-11T10:56:22+00:00">
                 2010/03/11 10:56
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               还没有学习animation包。求精确值用的代码是
               <code>
                integrate(f,0,1)
               </code>
               。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给邓一硕" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1250#respond" onclick='return addComment.moveForm( "div-comment-1250", "1250", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
            <ol class="children">
             <li class="comment byuser comment-author-yihui even depth-3" id="comment-1252">
              <article class="comment-body" id="div-comment-1252">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://yihui.name" rel="external nofollow">
                   谢益辉
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1252">
                  <time datetime="2010-03-11T11:30:09+00:00">
                   2010/03/11 11:30
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 索得斯呢……那么其实这个精确值也是近似值了，只不过近似方法不一样而已了。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-dengyishuo bypostauthor odd alt depth-3" id="comment-1253">
              <article class="comment-body" id="div-comment-1253">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/7969a4509f4093e8637e6bf70f3f7fc1?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  邓一硕
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1253">
                  <time datetime="2010-03-11T11:39:01+00:00">
                   2010/03/11 11:39
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 <code>
                  pnorm(1)-pnorm(0)
                 </code>
                 和
                 <code>
                  integrate(f,0,1)
                 </code>
                 结果一样都是0.3413447，都不是精确值。想更精确的话就用级数展开近似好了。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-yihui even depth-3" id="comment-1257">
              <article class="comment-body" id="div-comment-1257">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://yihui.name" rel="external nofollow">
                   谢益辉
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1257">
                  <time datetime="2010-03-11T12:54:51+00:00">
                   2010/03/11 12:54
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 倒也用不着这么“抠”啦，反正只要涉及计算机的东西几乎全都是近似……
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
            </ol>
            <!-- .children -->
           </li>
           <!-- #comment-## -->
           <li class="comment odd alt depth-2 parent" id="comment-1255">
            <article class="comment-body" id="div-comment-1255">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/edea739fb1cea19c3e57833e287e1011?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                colinisstudent
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1255">
                <time datetime="2010-03-11T12:43:02+00:00">
                 2010/03/11 12:43
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               pnorm应该是用一个很特殊的函数erf算得，具体的算法我也不知道。但是把一维积分转化成pnorm,把高维积分转化为高维正态分布函数是一个不错的思路，R里面的mvtnorm包里有pmvnorm函数。pnorm和pmvnorm都比别的近似方法快很多。如果无法转化，可以用Gauss-Hermit Quadrature或者拉普拉斯展开的方法求得近似的解，但是这两种方法应该都无法应付太高维的积分问题，估计2，3维顶天了
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给colinisstudent" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1255#respond" onclick='return addComment.moveForm( "div-comment-1255", "1255", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
            <ol class="children">
             <li class="comment byuser comment-author-dengyishuo bypostauthor even depth-3" id="comment-1258">
              <article class="comment-body" id="div-comment-1258">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/7969a4509f4093e8637e6bf70f3f7fc1?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  邓一硕
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1258">
                  <time datetime="2010-03-11T13:18:16+00:00">
                   2010/03/11 13:18
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 你说的erf是误差函数（又叫高斯误差函数），是个非基本函数，是用级数算的。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
            </ol>
            <!-- .children -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-1254">
          <article class="comment-body" id="div-comment-1254">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/edea739fb1cea19c3e57833e287e1011?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              colinisstudent
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1254">
              <time datetime="2010-03-11T12:30:34+00:00">
               2010/03/11 12:30
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             高维下计算积分除了MC还没有别的太好的方法。一个例子就是EM算法，积分维数经常相当的高，没办法，现在都是用的MC-EM的做法，先抽，再叠代优化似然函数，再抽，再叠代优化似然函数。我叫这种方法——抽风(疯)流
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给colinisstudent" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1254#respond" onclick='return addComment.moveForm( "div-comment-1254", "1254", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-yihui even thread-odd thread-alt depth-1 parent" id="comment-1259">
          <article class="comment-body" id="div-comment-1259">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              <a class="url" href="http://yihui.name" rel="external nofollow">
               谢益辉
              </a>
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1259">
              <time datetime="2010-03-11T13:42:57+00:00">
               2010/03/11 13:42
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             我再提供一份蒙特卡洛方法的学习材料吧：
             <a href="http://www.ceremade.dauphine.fr/~xian/coursMC.pdf" rel="nofollow">
              http://www.ceremade.dauphine.fr/~xian/coursMC.pdf
             </a>
             可算是讲义，很详尽。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给谢益辉" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1259#respond" onclick='return addComment.moveForm( "div-comment-1259", "1259", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2" id="comment-1268">
            <article class="comment-body" id="div-comment-1268">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/5c85d270961fb44f3b10bd684f0c60ee?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                <a class="url" href="http://www.sciencenet.cn/u/jiaojiaojing/" rel="external nofollow">
                 焦静
                </a>
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1268">
                <time datetime="2010-03-12T17:27:05+00:00">
                 2010/03/12 17:27
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               大致看了看，非常好的一份材料，不仅研究了Monte Carlo,还研究了Bayesian，非常值得大家看哦，不过缺点就在于这份PDF没有目录，所以我为这个材料做了一份目录，希望可以帮助大家学习：
              </p>
              <p>
               Monte Carlo Statistical Methods
               <br/>
               by Christian P.Robert
              </p>
              <p>
               context
              </p>
              <p>
               1  introduction
               <br/>
               1.1 Statistical Models                                5
               <br/>
               1.2 Likelihood Methods                                10
               <br/>
               1.3 Bayesian Methods                                  21
               <br/>
               1.4 Deterministic Numerical Methods                   28
               <br/>
               1.5 simulation versus numerical analysis:
               <br/>
               when is it useful?                                31
              </p>
              <p>
               2  Random Variable Generation                            35
               <br/>
               2.1 Basic Methods                                     37
               <br/>
               2.2 Beyond uniform distributions                      51
              </p>
              <p>
               3  Monte Carlo Intergration                              76
               <br/>
               3.1 introduction                                      77
               <br/>
               3.2 Classical Monte Carlo Integration                 80
               <br/>
               3.3 importance sampling                               87
               <br/>
               3.4 Acceleration Methods                              98
              </p>
              <p>
               4  Markov Chains                                         108
               <br/>
               4.1 Basic Notions                                     110
               <br/>
               4.2 Irreducibility                                    115
               <br/>
               4.3 Transience/Recurrence                             123
               <br/>
               4.4 Invariant Measures                                126
               <br/>
               4.5 Ergodicity and stationarity                       130
               <br/>
               4.6 Limit Theorems                                    134
              </p>
              <p>
               5  Monte Carlo Optimization                              139
               <br/>
               5.1 Introduction                                      140
               <br/>
               5.2 Stochastic Exploration                            142
               <br/>
               5.3 Stochastic Approximation                          173
               <br/>
               5.3.3 MCEM                                        195
               <br/>
               6 The Metropolis-Hastings Algorithm
               <br/>
               6.1 Markov Chain Monte Carlo                          197
               <br/>
               6.2 The Metropolis-Hastings Algorithm                 199
               <br/>
               6.3 A Collection of Metropolis-Hastings Algorithms    204
               <br/>
               6.4 Extensions                                        217
              </p>
              <p>
               7 The Gibbs Sampler                                      231
               <br/>
               7.1 General Principles                                232
               <br/>
               7.1.5 Hierarchical models                         253
               <br/>
               7.2 Data Augmentation                                 255
               <br/>
               7.3 Improper Priors                                   271
              </p>
              <p>
               8 Diagnosing Convergence                                 278
               <br/>
               8.1 Stopping the Chain                                279
               <br/>
               8.2 Monitoring Stationarity Convergence               282
               <br/>
               8.3 Monitoring Average Convergence                    290
              </p>
              <p>
               9 Implementation in Missing Data Models                  317
               <br/>
               9.1 First examples                                    319
               <br/>
               9.2 Finite mixtures of distributions                  340
               <br/>
               9.3 Extensions                                        354
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给焦静" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1268#respond" onclick='return addComment.moveForm( "div-comment-1268", "1268", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-taiyun even thread-even depth-1 parent" id="comment-1267">
          <article class="comment-body" id="div-comment-1267">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/beb15ba775be6241a4e5fc9e476e36c0?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              <a class="url" href="http://taiyun.cos.name/" rel="external nofollow">
               魏太云
              </a>
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1267">
              <time datetime="2010-03-12T14:47:18+00:00">
               2010/03/12 14:47
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             还是期待量子计算机的横空出世吧，到时候可不惧怕这些大计算量问题了。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给魏太云" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1267#respond" onclick='return addComment.moveForm( "div-comment-1267", "1267", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2" id="comment-1270">
            <article class="comment-body" id="div-comment-1270">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/edea739fb1cea19c3e57833e287e1011?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                colinisstudent
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1270">
                <time datetime="2010-03-14T03:41:24+00:00">
                 2010/03/14 03:41
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               我倒是觉得并不是所有的积分问题都值得硬上的，一些问题能够被简化，使得积分问题在相对比较低维的情况下进行。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给colinisstudent" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1270#respond" onclick='return addComment.moveForm( "div-comment-1270", "1270", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1" id="comment-1392">
          <article class="comment-body" id="div-comment-1392">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/34e0143fb09d6f1a82876586307be50c?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              huang shuai
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1392">
              <time datetime="2010-04-25T15:55:29+00:00">
               2010/04/25 15:55
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             <a href="http://www.stat.columbia.edu/~cook/movabletype/mlm/RMHMC_MG_BC_SC_REV_08_04_10.pdf" rel="nofollow">
              http://www.stat.columbia.edu/~cook/movabletype/mlm/RMHMC_MG_BC_SC_REV_08_04_10.pdf
             </a>
            </p>
            <p>
             刚在GELMAN博客上看到这篇文章，非常兴奋，这个东东将是蒙特卡洛计算的一个里程碑！
            </p>
            <p>
             P.S. 看到“有两位经济学家用该方法解作过1900多维计算。” 请问可有文献啊。这个维数有点吓人，如果他们成功了，那简直太值得一学了。。。有文献就拿出来共享嘛。。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给huang shuai" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1392#respond" onclick='return addComment.moveForm( "div-comment-1392", "1392", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-1694">
          <article class="comment-body" id="div-comment-1694">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/0cc73b0ff339b002fd626a989e261ea2?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              Zhangtao
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1694">
              <time datetime="2010-10-25T17:53:00+00:00">
               2010/10/25 17:53
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             非常喜欢MCMC，就是不太懂，希望大侠们多提供一些基础性的材料和讲义。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给Zhangtao" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1694#respond" onclick='return addComment.moveForm( "div-comment-1694", "1694", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1" id="comment-1723">
          <article class="comment-body" id="div-comment-1723">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/4f04a75519f1d0aee8d4af28bdd3ee3c?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              吴峰
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1723">
              <time datetime="2010-11-21T00:48:24+00:00">
               2010/11/21 00:48
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             受益了~下面的讨论也很精彩
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给吴峰" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1723#respond" onclick='return addComment.moveForm( "div-comment-1723", "1723", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1 parent" id="comment-1770">
          <article class="comment-body" id="div-comment-1770">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/5b8cd887d46c3894e895fbea7a46a3bf?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              W.C Zhu
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1770">
              <time datetime="2010-12-02T01:24:40+00:00">
               2010/12/02 01:24
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             如果能加点图例就更好了，像hit and miss方法，看图一目了然。再就是可以考虑再加上Importance sampling, MCMC等方法，感觉两种方法还是略显单调。最后说下版面问题：1.最上面因为页面右面有标签看起来还挺舒服，但是往下走走，感觉右边空白太大了（莫非是为了非宽屏电脑考虑？） 再就是留言评论呈阶梯状排列使得后面发言的人评论区很小，像谢老大写的代码被挤到一个小狭长的空间里了，都要托滚动条才能看完整，不知到能否改进（能不能让右边距也右移呢？）
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给W.C Zhu" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1770#respond" onclick='return addComment.moveForm( "div-comment-1770", "1770", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment byuser comment-author-yihui even depth-2" id="comment-1771">
            <article class="comment-body" id="div-comment-1771">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                <a class="url" href="http://yihui.name" rel="external nofollow">
                 谢益辉
                </a>
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-1771">
                <time datetime="2010-12-02T05:45:54+00:00">
                 2010/12/02 05:45
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               1、这个建议倒是提醒了我如何充分利用页面空间。往下看评论的时候页面右边确实显得有点空，此时不妨放一个浮动层，里面写点重要通知之类的，既能引起人的注意，又能改善空间布局。
               <br/>
               2、目前没办法，这种情况应该也不会经常出现。你可以双击代码，然后复制到别的地方看。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给谢益辉" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=1771#respond" onclick='return addComment.moveForm( "div-comment-1771", "1771", "respond", "2010" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="pingback odd alt thread-odd thread-alt depth-1" id="comment-2268">
          <div class="comment-body">
           Pingback：
           <a class="url" href="http://www.gaotao.name/cn/2010/03/improving-integrals-accuracy-by-m/" rel="external nofollow">
            提高Monte Carlo模拟定积分计算的精度 · 境由心生
           </a>
          </div>
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1" id="comment-3597">
          <article class="comment-body" id="div-comment-3597">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/586703045de709fe8330efea0b1494ac?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              小贝叶斯
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-3597">
              <time datetime="2013-01-18T14:55:55+00:00">
               2013/01/18 14:55
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             任何方法，在处理高维问题是，都会显得捉襟见肘。但是相比普通的数值积分，蒙特卡罗积分的收敛速度和维度无关。这是因为数值积分相当于在积分局域上均匀抽样，而蒙特卡罗积分则在“重要的区域”抽取多的样本，而不重要的区域抽取少的样本，这样的比较下，高下立判。当然了，现在的主要问题是如何从高维分布里抽样，这才是真正的难题！
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给小贝叶斯" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=3597#respond" onclick='return addComment.moveForm( "div-comment-3597", "3597", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-odd thread-alt depth-1" id="comment-4167">
          <article class="comment-body" id="div-comment-4167">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/5f90cd1a45ece3078680bf6dd16fd3b1?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              Crixus
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-4167">
              <time datetime="2013-05-06T20:20:51+00:00">
               2013/05/06 20:20
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             是不是E(f(X)) 应该= 1/n * ∫10f(x)dx=J???
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给Crixus" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=4167#respond" onclick='return addComment.moveForm( "div-comment-4167", "4167", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1" id="comment-4187">
          <article class="comment-body" id="div-comment-4187">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/6a65521e46d9ab80fd01be4e77284043?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              jessica_chu
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-4187">
              <time datetime="2013-05-11T21:19:28+00:00">
               2013/05/11 21:19
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             大哥，【进一步如果在区间[a,b]上有c≤g(x)≤d ，令f(y)=1d−cg(x)−c=1d−cg[a+(b−a)y]−c，则0≤f(y)≤1。此时，可以得到J′=∫bag(x)dx=S0J+c(b−a)。】这里少了括号好吗！
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给jessica_chu" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=4187#respond" onclick='return addComment.moveForm( "div-comment-4187", "4187", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="pingback odd alt thread-odd thread-alt depth-1" id="comment-5143">
          <div class="comment-body">
           Pingback：
           <a class="url" href="http://techtrace.org/?p=1399" rel="external nofollow">
            蒙特卡洛方法与定积分计算
           </a>
          </div>
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1" id="comment-5689">
          <article class="comment-body" id="div-comment-5689">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/b7d75199469978f21634a2fbc30fe976?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              圈圈
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-5689">
              <time datetime="2014-05-11T16:44:13+00:00">
               2014/05/11 16:44
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             不是有metropolis的重要性抽样方法吗，高维就不会稀疏了
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给圈圈" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=5689#respond" onclick='return addComment.moveForm( "div-comment-5689", "5689", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-odd thread-alt depth-1" id="comment-5694">
          <article class="comment-body" id="div-comment-5694">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/6b78af308888cb4e65be42c368bdf020?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              <a class="url" href="http://www.zhouleyu.com/life" rel="external nofollow">
               zhouleyu.com
              </a>
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#comment-5694">
              <time datetime="2014-05-12T14:38:07+00:00">
               2014/05/12 14:38
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             不是有metropolis的重要性抽样方法吗，高维就不会稀疏了
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给zhouleyu.com" class="comment-reply-link" href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/?replytocom=5694#respond" onclick='return addComment.moveForm( "div-comment-5694", "5694", "respond", "2010" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
        </ol>
        <!-- .comment-list -->
        <div class="comment-respond" id="respond">
         <h3 class="comment-reply-title" id="reply-title">
          发表评论
          <small>
           <a href="http://cos.name/2010/03/monte-carlo-method-to-compute-integration/#respond" id="cancel-comment-reply-link" rel="nofollow" style="display:none;">
            取消回复
           </a>
          </small>
         </h3>
         <form action="http://cos.name/wp-comments-post.php" class="comment-form" id="commentform" method="post" novalidate="">
          <p class="comment-notes">
           <span id="email-notes">
            电子邮件地址不会被公开。
           </span>
           必填项已用
           <span class="required">
            *
           </span>
           标注
          </p>
          <p class="comment-form-comment">
           <label for="comment">
            评论
           </label>
           <textarea aria-required="true" cols="45" id="comment" maxlength="65525" name="comment" required="required" rows="8">
           </textarea>
          </p>
          <p class="comment-form-author">
           <label for="author">
            姓名
            <span class="required">
             *
            </span>
           </label>
          </p>
          <p class="comment-form-email">
           <label for="email">
            电子邮件
            <span class="required">
             *
            </span>
           </label>
          </p>
          <p class="comment-form-url">
           <label for="url">
            站点
           </label>
          </p>
          <p class="form-submit">
          </p>
          <p style="display: none;">
          </p>
          <p style="display: none;">
          </p>
         </form>
        </div>
        <!-- #respond -->
       </div>
       <!-- #comments -->
      </div>
      <!-- #content -->
     </div>
     <!-- #primary -->
     <div class="sidebar-container" id="tertiary" role="complementary">
      <div class="sidebar-inner">
       <div class="widget-area">
        <aside class="widget widget_text" id="text-3">
         <h3 class="widget-title">
          关注统计之都
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            新浪微博
            <a href="http://weibo.com/cosname">
             @统计之都
            </a>
           </li>
           <li>
            人人网
            <a href="http://renren.com/cosname">
             @统计之都
            </a>
           </li>
           <li>
            Twitter
            <a href="http://twitter.com/cos_name">
             @cos_name
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-6">
         <h3 class="widget-title">
          微信公众平台
         </h3>
         <div class="textwidget">
          <img src="http://cos.name/wp-content/uploads/2013/04/qrcode-8cm.jpg"/>
          <p style="font-size:12px;margin-left:15px">
           微信号 CapStat
          </p>
          <p>
           我们将第一时间向您推送主站和论坛的精彩内容，以及统计之都的线下活动、竞赛、培训和会议信息。
          </p>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-8">
         <h3 class="widget-title">
          站内导航
         </h3>
         <div class="textwidget">
          <ul>
           <li class="page_item">
            <a href="http://cos.name/cn/">
             中文论坛
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/videos/">
             视频
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/salon/">
             统计沙龙
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/books/">
             图书出版
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/training/">
             教育培训
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/about/">
             关于我们
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/donate/">
             捐赠
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/chinar/">
             R语言会议
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_links" id="linkcat-2">
         <h3 class="widget-title">
          友情链接
         </h3>
         <ul class="xoxo blogroll">
          <li>
           <a href="http://stat.ruc.edu.cn" target="_blank" title="中国人民大学统计学院网站">
            中国人民大学统计学院
           </a>
          </li>
          <li>
           <a href="http://rucdmc.net">
            中国人民大学数据挖掘中心
           </a>
          </li>
          <li>
           <a href="http://birc.gsm.pku.edu.cn/" target="_blank">
            北京大学商务智能研究中心
           </a>
          </li>
          <li>
           <a href="http://sam.cufe.edu.cn/" target="_blank" title="中央财经大学统计与数学学院网站">
            中央财经大学统计与数学学院
           </a>
          </li>
          <li>
           <a href="http://tjx.cueb.edu.cn/" target="_blank" title="首都经济贸易大学统计学院网站">
            首经贸统计学院
           </a>
          </li>
          <li>
           <a href="http://www.shookr.com/">
            数客网大数据社区
           </a>
          </li>
          <li>
           <a href="http://www.xueqing.tv/" target="_blank" title="数据科学在线学习平台">
            雪晴数据网
           </a>
          </li>
          <li>
           <a href="http://iera.name/" target="_blank" title="IERA是一个旨在普及、传播和增进工业工程知识的非营利性网站">
            IERA（直通IE）
           </a>
          </li>
         </ul>
        </aside>
        <aside class="widget widget_categories" id="categories-2">
         <h3 class="widget-title">
          全部分类
         </h3>
         <label class="screen-reader-text" for="cat">
          全部分类
         </label>
         <select class="postform" id="cat" name="cat">
          <option value="-1">
           选择分类目录
          </option>
          <option class="level-0" value="925">
           cos访谈  (4)
          </option>
          <option class="level-0" value="659">
           图书出版  (5)
          </option>
          <option class="level-0" value="379">
           数学方法  (14)
          </option>
          <option class="level-1" value="381">
           分析与代数  (1)
          </option>
          <option class="level-1" value="380">
           概率论  (9)
          </option>
          <option class="level-1" value="382">
           随机过程  (5)
          </option>
          <option class="level-0" value="210">
           数据分析  (81)
          </option>
          <option class="level-1" value="203">
           多元统计  (3)
          </option>
          <option class="level-1" value="42">
           数据挖掘与机器学习  (42)
          </option>
          <option class="level-1" value="36">
           生物与医学统计  (17)
          </option>
          <option class="level-1" value="35">
           计量经济学  (4)
          </option>
          <option class="level-1" value="296">
           金融统计  (3)
          </option>
          <option class="level-1" value="303">
           风险精算  (7)
          </option>
          <option class="level-0" value="177">
           模型专题  (15)
          </option>
          <option class="level-1" value="38">
           回归分析  (10)
          </option>
          <option class="level-1" value="41">
           时间序列  (2)
          </option>
          <option class="level-0" value="784">
           每周精选  (24)
          </option>
          <option class="level-1" value="183">
           可视化  (9)
          </option>
          <option class="level-0" value="967">
           沙龙纪要  (3)
          </option>
          <option class="level-0" value="18">
           经典理论  (46)
          </option>
          <option class="level-1" value="37">
           抽样调查  (3)
          </option>
          <option class="level-1" value="4">
           统计推断  (26)
          </option>
          <option class="level-1" value="236">
           试验设计  (7)
          </option>
          <option class="level-1" value="39">
           非参数统计  (3)
          </option>
          <option class="level-0" value="1">
           统计之都  (279)
          </option>
          <option class="level-1" value="884">
           中国R会议  (2)
          </option>
          <option class="level-1" value="885">
           中国R语言会议  (5)
          </option>
          <option class="level-1" value="446">
           出国留学  (3)
          </option>
          <option class="level-1" value="179">
           推荐文章  (90)
          </option>
          <option class="level-1" value="3">
           新闻通知  (75)
          </option>
          <option class="level-1" value="263">
           统计刊物  (10)
          </option>
          <option class="level-1" value="174">
           网站导读  (40)
          </option>
          <option class="level-1" value="204">
           职业事业  (51)
          </option>
          <option class="level-1" value="213">
           高校课堂  (9)
          </option>
          <option class="level-0" value="178">
           统计计算  (28)
          </option>
          <option class="level-1" value="40">
           优化与模拟  (15)
          </option>
          <option class="level-1" value="43">
           贝叶斯方法  (6)
          </option>
          <option class="level-0" value="378">
           软件应用  (116)
          </option>
          <option class="level-1" value="44">
           统计图形  (36)
          </option>
          <option class="level-1" value="110">
           统计软件  (83)
          </option>
         </select>
        </aside>
        <aside class="widget widget_recent_comments" id="recent-comments-2">
         <h3 class="widget-title">
          最新评论
         </h3>
         <ul id="recentcomments">
          <li class="recentcomments">
           <span class="comment-author-link">
            fineboom
           </span>
           发表在《
           <a href="http://cos.name/2016/06/use-shiny-fleetly-set-up-visual-prototype-system/#comment-7317">
            利用shiny包快速搭建可视化原型系统
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            胡家新
           </span>
           发表在《
           <a href="http://cos.name/2016/06/r%e8%af%ad%e5%8d%83%e5%af%bb%e7%ac%ac%e4%b8%89%e6%9c%9f%ef%bc%9a%e5%bc%a0%e6%97%a0%e5%bf%8c%e7%a9%b6%e7%ab%9f%e7%88%b1%e8%b0%81%ef%bc%9f/#comment-7316">
            R语千寻第三期：张无忌究竟爱谁？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            <a class="url" href="http://www.zijiacha.com/category.php?id=6" rel="external nofollow">
             南糯山普洱茶
            </a>
           </span>
           发表在《
           <a href="http://cos.name/2016/06/r%e8%af%ad%e5%8d%83%e5%af%bb%e7%ac%ac%e4%b8%89%e6%9c%9f%ef%bc%9a%e5%bc%a0%e6%97%a0%e5%bf%8c%e7%a9%b6%e7%ab%9f%e7%88%b1%e8%b0%81%ef%bc%9f/#comment-7315">
            R语千寻第三期：张无忌究竟爱谁？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            J
           </span>
           发表在《
           <a href="http://cos.name/2016/05/value-of-the-reputation-from-the-data/#comment-7314">
            数据告诉你：高信誉的卖家应该收高价，还是收低价？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            <a class="url" href="http://gg" rel="external nofollow">
             Hilda
            </a>
           </span>
           发表在《
           <a href="http://cos.name/2013/01/drawing-map-in-r-era/#comment-7311">
            R时代，你要怎样画地图？
           </a>
           》
          </li>
         </ul>
        </aside>
        <aside class="widget widget_rss" id="rss-282869971">
         <h3 class="widget-title">
          <a class="rsswidget" href="http://cos.name/cn/topics/feed/">
           <img src="http://cos.name/wp-includes/images/rss.png"/>
          </a>
          <a class="rsswidget" href="http://cos.name/cn/topics/feed/">
           中文论坛新帖
          </a>
         </h3>
         <ul>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/417372/">
            处理时间数据和产生时间序列的问题
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/6790/">
            《统计陷阱》下载 （How to lie with statistics）
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/818/">
            统计学的世界（第五版）
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/16574/">
            class(x) 返回值值是AsIs,AsIs代表什么，有什么用处？
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/417366/">
            如何用R绘制一个分类算法的决策规则
           </a>
          </li>
         </ul>
        </aside>
        <aside class="widget widget_text" id="text-2">
         <h3 class="widget-title">
          登录/RSS
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <a href="http://cos.name/wp-admin/">
             登录
            </a>
           </li>
           <li>
            <a href="http://cos.name/feed/" title="使用 RSS 2.0 同步站点内容">
             文章
             <abbr title="Really Simple Syndication">
              RSS
             </abbr>
            </a>
           </li>
           <li>
            <a href="http://cos.name/comments/feed/" title="RSS 上的最近评论">
             评论
             <abbr title="Really Simple Syndication">
              RSS
             </abbr>
            </a>
           </li>
          </ul>
         </div>
        </aside>
       </div>
       <!-- .widget-area -->
      </div>
      <!-- .sidebar-inner -->
     </div>
     <!-- #tertiary -->
    </div>
    <!-- #main -->
    <footer class="site-footer" id="colophon" role="contentinfo">
     <div class="site-info">
      版权所有 © 2014 统计之都 | 由
      <a href="http://wordpress.org/">
       WordPress
      </a>
      构建 | 主题修改自
      <a href="http://wordpress.org/themes/twentythirteen">
       Twenty Thirteen
      </a>
     </div>
     <!-- .site-info -->
    </footer>
    <!-- #colophon -->
   </div>
   <!-- #page -->
   <p style="margin:0;padding:0;height:1px;overflow:hidden;">
    <a href="http://www.wumii.com/widget/relatedItems" style="border:0;">
     <img src="http://static.wumii.cn/images/pixel.png"/>
    </a>
   </p>
  </div>
 </body>
</html>